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Executive Summary 

When two or more Doppler weather radar systems are monitoring the same region, the Doppler 
velocities can be combined to form a three-dimensional (3-D) wind vector field thus providing for a more 
intuitive analysis of the wind field. A real-time display of the 3-D winds can assist forecasters in 
predicting the onset of convection and severe weather. The data can also be used to initialize local 
numerical weather prediction models. Two operational Doppler Radar systems are in the vicinity of 
Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS); these systems are 
operated by the 45th Space Wing (45 SW) and the National Weather Service Melbourne, Fla. (NWS 
MLB). Dual-Doppler applications were considered by the 45 SW in choosing the site for the new radar. 
Accordingly, the 45th Weather Squadron (45 WS), NWS MLB and the National Aeronautics and Space 
Administration tasked the Applied Meteorology Unit (AMU) to investigate the feasibility of establishing 
dual-Doppler capability using the two existing systems. This study investigated technical, hardware, and 
software requirements necessary to enable the establishment of a dual-Doppler capability. 

Review of the available literature pertaining to the dual-Doppler technique and consultation with 
experts revealed that the physical locations and resulting beam crossing angles of the 45 SW and NWS 
MLB radars make them ideally suited for a dual-Doppler capability. The dual-Doppler equations were 
derived to facilitate complete understanding of dual-Doppler synthesis; to determine the technical 
information requirements; and to determine the components of wind velocity from the equation of 
continuity and radial velocity data collected by the two Doppler radars. Analysis confirmed the suitability 
of the existing systems to provide the desired capability. In addition, it is possible that both 45 SW radar 
data and Terminal Doppler Weather Radar data from Orlando International Airport could be used to 
alleviate any radar geometry issues at the NWS MLB radar, such as the “cone of silence” or beam 
blockage. In the event of a radar outage at one of the sites, the multi-radar algorithms would provide 
continuing coverage of the area through use of the data from the remaining operational radar sites. 

There are several options to collect, edit, synthesize and display dual-Doppler data sets. These options 
include commercial packages available for purchase and a variety of freeware packages available from 
the National Center for Atmospheric Research (NCAR) for processing raw radar data. However, 
evaluation of the freeware packages revealed that they do not have sufficient documentation and 
configuration control to be certified for 45 SW use. Additionally, a T1 data line must be installed/leased 
from the NWS MLB office and CCAFS to enable the receipt of NWS MLB raw radar data to use in the 
dual-Doppler synthesis. Integration of the T1 data line into the Eastern Range infrastructure that will meet 
the security requirements necessary for 45 SW use is time-consuming and costly. 

Overall evaluation indicates that establishment of the dual-Doppler capability using the existing 
operational radar systems is desirable and feasible with no technical concerns. Installation of such a 
system represents a significant enhancement to forecasting capabilities at the 45 WS and at NWS MLB. 
However, data security and cost considerations must be evaluated in light of current budgetary 
constraints. In any case, gaining the dual-Doppler capability will provide opportunities for better 
visualization of the wind field and better forecasting of the onset of convection and severe weather events 
to support space launch operations at KSC and CCAFS. 
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1. Introduction 

When two or more Doppler weather radar systems are monitoring the same region, the Doppler 
velocities can be combined to form a three-dimensional (3-D) wind vector field. Such a wind field allows 
a more intuitive analysis of the 3-D winds, especially for users with little or no experience in deciphering 
Doppler velocities (Bousquet 2004). A real-time display of the wind field can assist forecasters in 
predicting the onset of convection and severe weather. The data can also be used to initialize local 
numerical weather models. Two Doppler radars are in the vicinity of Kennedy Space Center (KSC) and 
Cape Canaveral Air Force Station (CCAFS): the 45th Space Wing (45 SW) RadTec 43/250 radar and the 
National Weather Service (NWS) office in Melbourne, Fla. ( MLB) Weather Surveillance Radar - 1988 
Doppler (WSR-88D) radar. The orientation of these radars with KSC and CCAFS are shown in Figure 1. 
The 45th Weather Squadron (45 WS), NWS MLB and National Aeronautics and Space Administration 
(NASA) customers tasked the Applied Meteorology Unit (AMU) to investigate the feasibility of 
establishing dual-Doppler capability using these two systems. This task summarizes a literature review 
and consultation with experts to determine technical, hardware, and software requirements necessary to 
create a dual-Doppler capability. The AMU also investigated cost considerations and viable alternatives. 

Developing a dual-Doppler capability can provide a real-time horizontal wind field display to assist 
45 WS and NWS forecasters in predicting the onset of convection and severe weather events. These data 
can also be used to initialize a local mesoscale numerical weather prediction model. In addition, data 
combined from two or more radars mitigate radar geometry problems such as the cone of silence and 
beam blockage. (Lakshmanan et al. 2007) 



Figure 1. Orientation of the 45 SW radar and the NWS MLB radars 
with respect to KSC and CCAFS. 
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2. Principles of Dual-Doppler Radar Measurements 

With traditional single Doppler radar, only the radial component of the wind field can be measured. 
When two or more Doppler radars are surveying an area, the Doppler velocities can be combined to 
estimate the 3-D wind field. Each individual radar system scans according to its own particular strategy 
and time, therefore the velocities from each radar must be merged to a common coordinate system and 
reference time before wind synthesis can occur (Friedrich and Hagen 2004). 

The main objective of multiple-Doppler analyses, which combine observations collected by several 
neighboring Doppler radars to construct the 3-D wind field, is to ease the interpretation of airflow and to 
provide detailed information on the kinematics of observed storms that even a non-radar specialist can 
easily understand. Knowledge of storm kinematic structure can provide interesting insights on the storm 
dynamics through revealing mesoscale wind features that would be difficult to ascertain from just the 
radial velocity field. Another unique benefit of dual- or multiple-Doppler wind synthesis is to provide 
information on the structure of the vertical motion field (Bousquet, et al. 2008). 

2.1 Dual-Doppler Equation Derivations 

The AMU derived the dual-Doppler equations to facilitate complete understanding of dual-Doppler 
synthesis, to determine the technical information requirements, and to determine how to find the three 
components of wind velocity from the continuity equation and radial velocity data from two Doppler 
radars (Armijo 1969). The derivations are available in the Appendix of this report for the interested 
reader. The dual-Doppler equations were compiled from a combination of the methods of Armijo (1969), 
Obrien (1970), Lhermitte and Miller (1970) and Carey (2005). The geometry of the radial distance 
between two radars and a point in space is shown in Figure 2. The typical beam geometry for a single 
radar is shown in Figure 3. 



Figure 2. The geometry of the radial distance (R) between radars and a point in space. The 
target point of the radars is at P(x, y, z). Radar 1 is at point (xi, yi, z,) and radar 2 is at point 
(x 2 , y 2 , z 2 ). Adapted from Figure 9.3 in Doviak and Zrnic (1993). 
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target 



Figure 3. The typical beam geometry for a single radar. Adapted from 
Figure 6.2 in Rinehart (2004). 

2.2 Dual-Doppler Area Geometry 

When planning a dual-Doppler capability, the specific area of interest must be studied to determine if 
the locations of the radars will allow an accurate wind field calculation. Dual-Doppler coverage depends 
on 

• The minimum spatial resolution needed to resolve the weather phenomena of interest, 

• The largest acceptable error in horizontal velocity, and 

• The distance between the radars (Davies-Jones 1979). 

• The orientation of the radars relative to the area for which dual-Doppler coverage is desired. 


2.2.1 Optimal vs. Actual Arrangement of Radars 

To begin the analysis, unit vectors of the two radial velocities intersect an angle, p (Davies-Jones 
1979, Friedrich and Hagen 2004). Using the schematic in Figure 2 as an example, p is the angle between 
Ri and R 2 at P. The error variances, a £ and <r^, of u- and v- wind components are related to the Doppler 
mean velocity error variances of the two individual radars, of and a | , by (Davies-Jones 1979): 


<r + CT~ 

U __ V 

_2 2 

<T, + G -y 


= CSC 2 /?. 


The optimal P that will provide the most accurate wind synthesis is 90° (Beck 2004). Experience at 
the National Severe Storms Laboratory (NSSL) has shown that a P of > 30° and < 150° is adequate for 
accurate wind synthesis (Davies-Jones 1979). Outside of this range, a large component of both beams is 
oriented in the same direction. This will cause large errors in the velocity components derived from the 
dual-Doppler wind synthesis. The calculated beam crossing angle between the 45 SW and the NWS MLB 
radars for launch complexes 17A (Figure 4a) and 41 (Figure 4b) on CCAFS are 67.9° and 52.2°, 
respectively, both within the 30° to 150° range recommended by Davies-Jones (1979). 
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Figure 4. a) The beam crossing angle between the 45 SW and NWS MLB radars for launch complexes 
17A (67.9°) and b) 41 (52.2°) on CCAFS. The yellow pins show the locations of the radars. The red lines 
labeled R| and R 2 are analogous to the similarly labeled lines in Figure 2. 


The spatial resolution of the dual-Doppler coverage area is related to the baseline, or horizontal 
distance between two radars. Davies-Jones (1979) and Friedrich and Hagen (2004) proposed that the 
optimal distance between two radars would range between 23.2 NM (43 km) and 43.7 NM (81 km). 
Examination of Figure 3 shows that the baseline between the 45 SW and NWS MLB radars is 23 NM 
(42.6 km), at the low end of the optimal separation distance. The beam crossing angles over KSC and 
CCAFS and the baseline of the 45 SW and NWS MLB radars makes them ideally suited for a dual- 
Doppler capability. A critical element of site location is the technical trade-off between velocity 
component accuracy and spatial resolution (Davies-Jones, 1979). This placement was considered by the 
45 WS during the original 45 SW radar site location determination. 

2.2.2 Dual-Doppler Area for 45 SW and NWS-MLB Radars 

The next step in the analysis is to calculate the total coverage area of the dual-Doppler analysis. 
Davies-Jones (1979) and Friedrich and Hagen (2004) show that the total coverage area is the intersection 
of the area (A|) defined by the upper limits on velocity error variance (Figure 5a) and the area (A 2 ) 
defined by the maximum range of the radars (Figure 5b). The area A t is made up of two circular areas 
called dual-Doppler lobes. In the formula for Ai (Davies-Jones 1979) shown in Figure 5a, d is half of the 
baseline (23 NM/2 =11.5 NM) and p is the minimum beam crossing angle in radians (71/6 rad or 30°). In 
the formula for A 2 in Figure 5b, R is the maximum range from any point in area A, to either of the radars. 
R for the 45 SW and NWS MLB radar pair is given by (Davies-Jones 1979) 

2*d*cscp = 2*( 11.51 NM)*csc(30°) = 46. 1 NM (85.3 km). 

Area A 2 is the maximum area that can be covered by two radars at an elevation angle of 0°. With 
increasing elevation, area A 2 decreases (Friedrich and Hagen 2004, Davies-Jones 1979). The dual- 
Doppler coverage area is defined as the area common to both A, and A 2 (Davies-Jones 1979): 

A n (/3,R)=A x {j})nA 2 (R) 
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Figure 5. a) Area A| that shows the dual-Doppler lobes as stippled areas bounded by P = 30° (;t/6 rad) 
(Figure la in Davies-Jones 1979), and b) Area A 2 , shown as a stippled area that lies within distance R of 
both radars (Figure lb in Davies-Jones 1979). 

In a dual-Doppler analysis, the resolution of the resulting merged data is defined as s = R7iA/180, 
where A is the half power beam width in degrees of the radars (Davies-Jones 1979). The technical 
specifications for the NWS MLB WSR-88D and the 45 SW Radtec 43/250 radars are shown in Table 1 . 
Using these values, the resolution of the 45 SW and NWS MLB radar pair is 

s = (46.1 NM)*ji*( 0.95')/180 = 0.76 NM (1.41 km). 

Using the formula for Ai shown in Figure 5a, the coverage area A! of the radar pair is 3,140.2 NM 2 
(10,770.7 km 2 ). Area A| for the 45 SW and NWS MLB radars is shown in Figure 6a. The maximum area 
that can be covered by the radar pair at an elevation angle of 0° is calculated using the formula for area A 2 
shown in Figure 5b. This area is shown in Figure 6b and is 4,565.6 NM 2 (15,659.7 km 2 ). The size of area 
A] fits into the size of area A 2 and is shown in Figure 7. Therefore the area common to both A, and A 2 for 
the 45 SW and NWS MLB radar pair is 3,140.2 NM 2 (10,770.7 km 2 ). However, the maximum area that 
can be covered by a radar pair decreases with increasing elevation angle (Davies-Jones 1979, Friedrich 
and Hagen 2004). 
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Table 1 . Technical specifications for the NWS MLB WSR-88D and 45 SW 
Radtec 43/250 radars. 

Parameter 

NWS-MLB WSR-88D 
(Salazar et al. 2009) 

45 SW Radtec 43/250 
(Short 2008) 

Wavelength 

10.5 cm 

5.32 cm 

Frequency 

S-band, 2.7 - 3.0 GHz 

C-band, 5.6-5.65 GHz 

Peak Power 

750 kW 

250 kW 

Antenna diameter 

8.5 m 

4.3 m 

Beam width 

0.95° 

0.95° 

Rotation Rate 

6 rpm 

6 rpm 

PRF 

318- 1304 s’ 1 

200 - 2000 s' 1 

Pulse Length 

1.57 and 4.7 ps 

0.4 to 16.7 ps 

Gain 

45.5 dB 

44 dB 

Latitude 

28.1131 N 

28.3938 N 

Longitude 

80.6541 W 

80.9510 W 



45 SW radar 


45 SW radar 


Figure 6. Same as Figure 4 except for the 45 SW and NWS MLB radars, where a) is area A, and 
b) is area A 2 . The yellow pins show the locations of the radars. The white circles in a) are the 
dual-Doppler lobes within which the angles subtended by the radials from the radars lie between 
30° and 150°. The intersecting area between the two circles along the radar baseline is not part of 
the dual-Doppler area because the beam crossing angles in this area are < 30° or > 1 50°. The 
white outline in b) shows area A 2 , the outline of points that lie within a distance of 46.1 NM (85.3 
km) of each radar. 




Figure 7. The intersection of areas A| and A 2 for the 45 SW and 
NWS MLB radars. The yellow pins show the locations of the 
radars. 
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Hardware and Software 

There are several options to collect, edit, synthesize and display dual-Doppler data sets. The 45 SW 
currently uses the Interactive Radar Information System (IRIS™) software package by Vaisala to display 
their radar data. The IRIS™ software has an add-on product, called NDOP, that provides capability for 
Dual-Doppler wind fields based on radial wind inputs from two Doppler radars that can ingest WSR-88D 
data. Ingesting WSR-88D data also provides the opportunity to eliminate the cone of silence, an important 
issue for thunderstorms approaching from the southwest (a frequent trajectory in the summer). The 
license includes the ability to make mosaics, or composites, of radar products from multiple sites. The list 
cost is $16,000. This is a one-time cost for an add-on license to the existing IRIS software used by the 45 
SW. An example of an IRIS™ NDOP product display is shown in Figure 8. Also a variety of freeware 
packages are available from the National Center for Atmospheric Research (NCAR) for processing raw 
radar data, but these packages do not have the thorough documentation and stringent configuration 
control needed in order to be certified for 45 SW use. 

Regardless of software choice, a T1 data line must be installed between the Morrell Operations 
Center (MOC) on CCAFS and NWS MLB to enable the receipt of NWS MLB raw radar data needed for 
the dual-Doppler synthesis. This line, with additional costs, could also be used to send the 45 SW radar 
data to NWS MLB. In this case, NWS MLB can use the multi-sensor, multi-radar processing options via 
the Warning Decision Support System Integrated Information (WDSS-II) system for viewing the dual- 
Doppler radar data and then pass this information back to the 45 SW. The AMU could then create a way 
to view the dual-Doppler data using Meteorological Interactive Data Display System (MIDDS), 
Advanced Weather Interactive Processing System (AWIPS), or other display techniques. Due to the 
proprietary nature of the vendor costs of installing and maintaining the T1 line, performing drawing and 
configuration changes, and assuring information security functions, only a lower bound on the cost could 
be obtained from the 45 SW. The cost and time estimate given was a minimum of $150,000 and 6 to 12 
months for this project. A monthly maintenance cost of the T1 line was estimated to be $ 1 000. 
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Figure 8. An example of an NDOP Display for strong northerly 
winter monsoon flow over China near Hong Kong, colliding with 
easterlies over the sea. The “+” points show regions where the 
minimum crossing angle between radars was met, but there were 
no weather targets, so a wind could not be calculated. The green 
and gray contours indicate ~300 m (-1000 ft) elevation intervals 
(Vaisala IRIS™ Product and Display Manual 2006). 
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3. Summary and Conclusions 

This study investigated technical, hardware, and software requirements necessary to enable the 
establishment of a dual-Doppler capability using the 45 SW Doppler and the NWS MLB WSR-88D 
radars. Use of these two Doppler radars to provide a 3-D display of the wind field will help forecasters to 
better predict the onset of convection and severe weather in support of space launch operations at KSC 
and CCAFS. The data can also be used to initialize a local mesoscale numerical weather prediction model 
to improve forecast output. 

3.1 Dual-Doppler Analysis Feasibility 

Evaluation of the physical locations and resulting beam crossing angles of the 45 SW and NWS MLB 
radars make them ideally suited for a dual-Doppler capability. The AMU derived the dual-Doppler 
equations to facilitate complete understanding of dual-Doppler synthesis, to determine the technical 
information requirements, and to determine how to find the three components of wind velocity from the 
equation of continuity and radial velocity data collected by two Doppler weather radars. 

There are several options to collect, edit, synthesize and display dual-Doppler data sets. In addition to 
the commercially available Vaisala NDOP package available for purchase, there are a variety of freeware 
packages available including WDSS-II software available via license from the University of Oklahoma 
and from NCAR for processing raw radar data. Unfortunately, the freeware packages do not have 
sufficient documentation and configuration control to be certified for 45 SW use. In any case, a T1 data 
line must be installed between the MOC and NWS MLB to enable the receipt of NWS MLB raw radar 
data to use in the dual-Doppler synthesis. The installation of a T1 data line that will meet the security 
requirements necessary for 45 SW use is time-consuming and costly. 

3.2 45 WS Interests 

Possible applications for utilization of dual-Doppler data by the 45 WS are as follows. Note that due 
to the distance of the NWS MLB radar from sites on KSC/CCAFS, the beam altitude of the lowest beam 
may be too high to see some of the low level features suggested below: 

• Easier visualization of wind fields for launch operation, pre-launch operations, wind 
warnings and wind forecasting. 

• Easier detection of threatening convective winds. 

• Better detection of convergence lines for forecasting convection initiation. 

• Detection of low level helicity for severe weather prediction. 

• Detection of winds near offshore for possible improved wind warnings. 

• Initiation of short forecast local numerical models, especially the wind field over the 
ocean that is currently not measured well. 

• Replacement of lost inland wind tower data to support 2 1 sl Century Spaceport or Air 
Force Smart Operations 21 s1 Century (AFSO-21) initiative. 

3.3 NWS MLB Interests 

NWS MLB agrees that the siting of the two radars is ideal for use of dual-Doppler techniques at 
KSC and CCAFS. Once the 45 SW and NWS MLB radar data co-exist. NWS MLB can use the multi- 
sensor, multi-radar processing options via the Warning Decision Support System Integrated Information 
(WDSS-II) system for viewing. The WDSS-II is the second generation of a system of tools for the 
analysis, diagnosis and visualization of remotely sensed weather data. WDSS-II has the capability to 
merge multiple-radar data into four-dimensional grids including Terminal Doppler Weather Radar 
(TDWR) data. It is possible that both 45 SW radar data and TDWR data from Orlando International 
Airport could be used to alleviate radar geometry issues at the NWS MLB radar, such as the cone of 
silence or beam blockage. In addition, in the event of a radar outage at one of the sites, the multi-radar 
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algorithms would provide continuing coverage of the area through use of the data from the remaining 
operational radar sites. 

3.4 Future Work 

Accurate near-term forecasting of convective activity and severe storms is of continuing interest to 
all customers. As has been described, dual-Doppler capability can provide significant improvement to that 
forecasting capability. To illustrate the level of improvement that may be achieved, the AMU suggests the 
dual-Doppler data be subjected to a proof-of-concept activity by which wind data is validated to provide 
baseline confidence. Specifically, dual-Doppler data can be compared to balloon and/or wind profile data 
to evaluate and quantify the observed accuracy. These proof-of -concept evaluations could be performed 
on a variety of “typical” days including a summer day and a winter day when no convection or frontal 
zones are present and the normal sea breeze location is known. Once the data is validated and the 
capability has been confirmed, merged scans from the existing radars during known periods of convective 
activity and storms, such as that observed Oct 9 -10, 2011 (Columbus Day Storm) can be analyzed to 
determine if the dual-Doppler capability would have aided in the near-term forecasting. 
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Appendix 

Derivation of the Dual-Doppler Equations 

The following derivations were compiled from Armijo (1969), Obrien (1970), Lhermitte and Miller 
(1970) and Carey (2005). Before the Doppler radar data can be used to perform dual-Doppler wind 
retrieval, they must be converted from polar to Cartesian coordinates as follows (Goh and Holt 2004). 

Referencing Figure 2, the distance between the target and radar 1 is: 


+{y-y i) 2 +( 2 -Zi) 2 

(A-l) 

The distance between the target and radar 2 is: 


R 2 =^{x-x 2 f +(y-y 2 f +{z-z 2 f 

(A-2) 

Squaring equations A-l and A-2 above and taking the derivative of the squared radial velocities with 
respect to time will allow the distances to be converted to velocities that can then be solved to find the 
horizontal wind components: 

R \ =(*-*i) 2 +(y-y\Y + { z ~ z \Y 

(A-3) 

R l =(*-* 2 ) 2 +(y-yiY + (z- z 2 ) 2 

(A-4) 

2R t = 2(x - x, )— + 2 (y - y ] )— + 2 (z - z, }^- 

1 dt V l 'dt u dt V ir dt 

(A-5) 

2 r 2 = 2(x - x 2 )— + 2(y - y 2 )— + 2(z - z 2 )— 

2 dt V 2J dt V 2, dt V 2> dt 

(A-6) 


Divide equation A-5 by 2R| and substitute Vi for , u for — , v for — , and w + Vj for — to get: 
M J dt dt dt dt 


V x = u 


(x-x.) (y-Ti) / ,,\{ z - z \) 

v 17 + v — — — + {w + V l )- 17 


R, 


R, 


R, 


(A-7) 


Divide equation A-6 by 2Ri and substitute V? for ^2 . , u for — , v for — , and w + V, for — to get: 

dt dt dt dt 


A P P V ' 7 p 

A ^ A -s I \ -) 


(A-8) 


Referencing Error! Reference source not found., write equations A-7 and A-8 in terms of the radars’ 
elevation (a) and azimuth (p) angles using the radar geometry as follows (Note: m denotes a generic radar 
number). 

~ * -■ = cosa m sin <t> m , ^ y - y — = cosa m cos^ m , and = sin a m . 

K„, K n, K„, 


Using the above substitutions, equation A-7 becomes: 

Fj =«cosa l sin^ 1 -i-vcosa, cos$ +(w + fj)sinor, . (A-9) 

Equation A-8 becomes: 

V 2 = ucosa 2 sin + vcosa 2 cos^ 2 +(w + V,)sma 2 . (A- 10) 
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Solve equation A-9 for v: 

v = seca, sec^,(-MCOsa, sin^, -(w+ F, )sina, + V ] ). 

Substitute equation A-l 1 into equation A- 10: 

V 2 =u cos a 2 sin <j> 2 + (sec a, sec fa (- u cos a, sin fa - (w + V, )sin a, + V ] )) 

* cos a 2 cos <f >2 ■+■ (w 4- V, )sin a 2 . 

Expand equation A- 12 and multiply both sides of equation A- 12 by cos Pi: 

V 2 cosfa = wcos« 2 sin cos$ -wcosa, sin$ seca, sec^, cosa 2 cos^ 2 cos^, 

- (w + V, )sina, seca, sec^, cosa 2 cos^ 2 cos$ + V x seca, sec^, cosa 2 cos$> cos$ 
+ (w + V, )sina 2 cos^, . 


Simplify equation A-l 3: 

F 2 cos$ = ucosa 2 sin cos$ - wsin^, cos^ 2 cosa 2 
- (w + V,)tana, cosa 2 cos^ 2 + V x seca, cosa 2 cos^ 2 
+ (w + V, )sina 2 cos^ 

Solve equation A- 14 for u: 

u cos a 2 (sin fa cos - sin <j> 2 cos <j ) ) = V x sec a, cos a 2 cos fa 

- V 2 cos^, - (w + Vjtana, cosa 2 cos^ 2 + (w + V,)sina 2 cos^, 

V t sec a, cos a 2 cos fa F 2 cos fa 

cos a 2 (sin fa cos fa - sin <j) 2 cos fa ) cos a 2 (sin fa cos fa - sin fa cos fa ) 


1 

(w + V, ) 

Itana, cos oc 2 cos^ 2 

i 

(w + V, ) 

lsina 2 cos^. 

cosa 2 ! 

(sin 0, cos - sin (j) 2 cos </> x ) 

i 

1 cosa 2 l 

(sin fa cos fa - sin fa cos fa ) 


1 

Vj COS <t> 2 v 2 COS fa 

(sin <j) x Co: 
(w + V t j 

5 fa - sin 02 COSfa ) 

cos tana, - co: 

cos a, 
>fa tana 2 

cosa 2 

(sin <j) x cos (j) 2 - sin ^ 2 cos ) 



(A-l 1) 


(A-l 2) 


(A-13) 


(A- 14) 


(A-l 5) 


(A- 16) 


(A- 17) 


Substitute the trigonometric identity sin(u - v) = sin u cos v - cos u sin v in the denominator of equation 
A- 17 (Lhermitte and Miller 1970): 


u - 


1 

sin (0i ~<f> 2 ) 


V l COS <t> 2 
cosa. 


F 2 cos fa 
cosa 2 


( ■J cos ^ tana, -cos^, tana 2 

' l ,JI -hto-^) 


sint 


(A- 18) 


Solve for v by repeating the above steps. Solve equation A- 10 for u: 

u = seca 2 csc$,(- vcosa 2 cos^ 2 -(w+ F,)sina 2 + ^)- 


(A- 19) 
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Substitute equation A- 19 into equation A-9: 

V. = csct/)J-vcnsa 2 cos<f> 2 - (w + Fj )sin a 2 +F 2 )]cosa! sin^, 


Fj = [seca 2 csc^ 2 (- vcosa 2 cos^ 2 -( 

-f vcoscifj cos^, +(w + V t )s\na { . 

and multiply both sides of equation A- 12 by sin p 2 - 


Expand equation A-20 and multiply both sic 

Fj sin ^ 2 = -v cos a 2 cos ^ 2 cos a l sin fa sec a 2 esc (p 2 sin (p 2 +v c 

(„r _l \r ~ ~ sin^, sin^ 2 + F 2 secar 2 csc^ 2 cos a, 


sec ar 2 esc (f> 2 sin (j) 2 + v cos a ] cos ^ sin fa 

* — -Vj sin^, sin ^ 2 


Fj sin ^ 2 = -v cos a 2 cos fa cos a x sin fa 
-(w + V t )sina 2 seca 2 csc^ 2 cos a, 

+ (w + V t )sin a ] sin fa 2 . 

Simplify equation A-21: 

Fj sin ^ 2 = -v cos fa 2 cos a, sin + v cos a ] cos ^ sin fa 2 

- (w + V t )tanar 2 cos a, sin^, + V 2 seca 2 cos a, sin fa 
(w + V t )sina, sin^ 2 - 


Solve equation A-22 for v: 

z, (cos fa 2 sin fa 


vc 


V : 


cos a. 


V 2 sec a 2 Co: 
:j (sin fa cos 


- sin fa 2 cos fa) = V 2 sec a 2 cos a x sin fa 
\ )tan a 2 cos a, sin fa + (w + V t )sin a, sin fa 2 

a j/ a 


- sin <f > 2 1 


1 

(w + V, ) 

ltana 2 cos a, sin$ 

1 

(w + V,) 

Isin sin /? 2 

cos a ! 1 

(sin <j> x cos (f> 2 - sin (f> 2 cos <j > x ) 

j 

1 cos a, 1 

(sin fa cos (f) 2 - sin (j) 2 cos fa ) 


i 

F 2 sin <j\ 

Fj sin (f) 2 

- sin <j> 2 cos <f >\ ) 

cos a 2 

cos a, 


( w + v j s ' n $ tan a 2 ~ s ' n ^2 tan a \ 

1 (sin <j> x cos <j> 2 - sin <j> 2 cos ) 


Substitute the trigonometric identity sin(u - v) = sin u cos v - cos u sin v in the denominator 
A-25 (Lhermitte and Miller 1970): 

_ 1 V 2 sin <f\ F,sin^ 2 

i($ -<t> 2 ) cos a 2 cosa. 


(A-20) 


(A-21) 


(A-2: 


(A 


sim 



Now there are two equations, A-18 and A-26, with four unknowns: u, v, w, and V,. Therefore, equations 
A- 18 and A-26 must be supplemented with the equation of continuity: 

dt 


d(pu) | d{pv) | d{pw) Q 


dx 


dy 


dz 


(A-27) 


Expand equation A-27 to get: 

dp du dp dv dp dw dp . 
h p h U h p h V h p 1- W — 0 

dt dx dx dy dy dz dz 


(A-28) 
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Assume density, p, is constant at a particular level, z, and simplify equation A-28: 

du dv dw dp . 
p — + p — + p — + w— = 0. 
dx dy dz dz 

Rearrange equation A-29: 


dw dp 

p 1 -w— - - 

dz dz 


du dv 

p — + p — 

dx dy 


(A-29) 


(A-30) 


„ . . d\pw) dw dp . . 

Substitute ^ = p — + w— into equation A-30: 

dz dz dz 


d(pw) _ 

dz 


du dv'' 

p — + p — 

dx dy 


(A-3 1 ) 


Now, equation A-3 1 can be integrated from one level to the next for each level: 

t h d(pw) e h ( du dv 

l*- 

/ _ 

Mb ~Ma=~ 



( du 

dv } 

Az - p 

( du 

dv) 


p 

1 

1 

Az 


^dx 

dy) 

b 

dx 

dy) 

a 


(A-32) 

(A-33) 


Equation A-33 can be solved for vertical velocity, w, at every height by first assuming w=0 at the surfc 
and then iterating until the mean difference between the horizontal and vertical divergence is less than 
some predetermined threshold. 

Finally, V t can be estimated by using standard V r Z relationships (Doviak et al. 1976): 

/ \0.4 

V t = 2.65Z 0114 — ms~' . 

\r 

The IRIS NDOP software estimates V, using the following relationships (IRIS 2006): 

Above the melting level (snow and graupel): V t = 0.8Z 006 [wj"' 

Below the melting level (rain): V t = 2.70Z 0 " [m -1 ] 


(A-34) 

(A-3 5) 
(A-3 6) 
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List of Acronyms 


45 SW 

45th Space Wing 

NCAR 

National Center for Atmospheric 

45 WS 

45th Weather Squadron 


Research 

AMU 

Applied Meteorology Unit 

NDOP 

IRIS product that calculates dual- 
Doppler wind fields based on 

AWIPS 

Advanced Weather Interactive 


radial wind inputs from two 


Processing System 


Doppler radars 

CCAFS 

Cape Canaveral Air Force Station 

NM 

Nautical Mile 

IRIS™ 

Interactive Radar Information 

NSSL 

National Severe Storms 


System 


Laboratory 

KSC 

Kennedy Space Center 

NWS 

National Weather Service 

MIDDS 

Meteorological Interactive Data 

TDWR 

Terminal Doppler Weather Radar 


Display System 

WDSS-II 

Warning Decision Support 

MLB 

Melbourne, Fla. - 3 letter 
identifier 


System Integrated Information 


WSR-88D 

Weather Surveillance Radar - 

MOC 

Morrell Operations Center 


1988 Doppler 

NASA 

National Aeronautics and Space 
Administration 
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NOTICE 


Mention of a copyrighted, trademarked or proprietary product, service, or document does not constitute 
endorsement thereof by the author, the AMU, the National Aeronautics and Space Administration, or the 
United States Government. Any such mention is solely for the purpose of fully informing the reader of 
the resources used to conduct the work reported herein. 
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